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Abstract. Noncommutativc rational functions appeared in many contexts 
in system theory and control, from the theory of finite automata and for- 
mal languages to robust control and LMIs. We survey the construction of 
noncommutative rational functions, their realization theory and some of their 
applications. We also develop a diff'erence-differential calculus as a tool for 
further analysis. 

1. Introduction 

Noncommutative rational functions first appeared in system theory in the con- 
text of recognizable formal power series in noncommuting indeterminates in the 
theory of formal languages and finite automata; see Kleene |36| . Schiitzenberger 
[48l |49] , and Fliess [23l [24l [25] (where the motivation comes also from applications 
to certain classes of nonlinear systems), and Berstel-Reutenauer [12] for a survey. 
In particular, noncommutative rational functions admit a good state space real- 
ization theory. More recently, state space realizations of rational expressions in 
Hilbert space operators (modelling structured possibly time varying uncertainty) 
have figured prominently in work on robust control of linear systems, see Beck [9], 
Beck-Doyle-Glover [10], Lu-Zhou-Doyle g^. 

Another important application comes from the area of Linear Matrix Inequali- 
ties (LMIs); see, e.g., Nesterov-Nemirovski [42], Nemirovski ^1], Skelton-Iwasaki- 
Grigoriadis |50j. As it turns out, most optimization problems appearing in sys- 
tems and control are dimension- independent, i.e., the natural variables are matri- 
ces, and the problem involves rational expressions in these matrix variables which 
have therefore the same form independent of matrix sizes; see Helton [29] . Helton- 
McCullough-Putinar-Vinnikov [30]. Realizations of rational functions in noncom- 
muting indeterminates are exactly what is needed here to convert (numerically 
unmanageable) rational matrix inequalities into (highly manageable) linear matrix 
inequalities, see Helton-McCullough-Vinnikov [31]. 

Last but not least, in many situations one can establish a commutative result 
by "lifting" to the noncommutative setting, applying the noncommutative theory, 
and then "descending" again to the commutative situation. Some examples are: 
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• The classical paper of Fornasini-Marchesini [27] establishing a state space 
realization theorem for rational functions of several commuting variables. 

• The results of Ball-Kaliuzhnyi-Verbovetskyi VP on conservative dilations of 
various classes of (commutative) multidimensional systems. 

• The theorem of Kaliuzhnyi-Verbovetskyi-Vinnikov |34| showing that the 
singularities of a matrix- valued rational function of several commuting vari- 
ables which is regular at zero coincide with the singularities of the resolvent 
in any of its Fornasini-Marchesini realizations with the minimal possible 
state space dimension. This implies, in particular, the absence of zero-pole 
cancellations in a minimal factorization. 

The goal of this paper is two-fold. First, we survey the basic concepts of the 
theory of noncommutative rational functions, and their realization theory. Second, 
we develop a diiference-differential calculus for noncommutative rational functions. 
This is a new powerful tool for the needs of system theory and beyond. It is a spe- 
cial instance of the general theory of noncommutative functions which are defined 
as functions on tuples of matrices of all sizes satisfying certain compatibility con- 
ditions as we vary the size of matrices (they respect direct sums and simultaneous 
similarities); see Kaliuzhnyi-Verbovetskyi-Vinnikov |35] . 

It is important to notice that the NCAlgebra software, 

http://www.math.ucsd.edu/~ncalg 

implements many symbolic algorithms in the noncommutative setting; see [30] for 
examples, guidance, and detailed references. 

2. Noncommutative rational functions 

We first formally introduce noncommutative rational functions; this involves 
some non-trivial details since unlike the commutative case, a noncommutative ra- 
tional function does not admit a canonical coprime fraction representation. We fol- 
low Kaliuzhnyi-Verbovetskyi-Vinnikov [34], to which we refer for both details and 
references to extensive algebraic literature; we only mention Amitsur [3], Bergmann 
Cohn [T5l[T7] for some of the original constructions, and Rowen (47j Chapter 
8], Cohn [16l[20] for good expositions. 

We start with noncommutative polynomials in d noncommuting indeterminates 
zi, . . . , Zd over a field K. E.g., a noncommutative polynomial of total degree 2 in 2 
indeterminates zi , Z2 is of the form 

p ^ a + (3zi + 7Z2 + Szf + eziZ2 + Cz2Zi + rjzl, 

where the coefficients a, /3, 7, 5, e^Q,rj ^ K. The general form of a noncommutative 
polynomial is 

Here J- a denotes the free semigroup with d generators (letters) gi, . . . , g^; elements 
of Fd are arbitrary words w = gi^ ■ ■ ■ gi^ and the semigroup operation is concatena- 
tion; the neutral element is the empty word 0, and \'w\ = ^ is the length of the word 
w. We use noncommutative multipowers = z^^ • • • Zi-^ . In the example above, 

P0 = a, Pg^ = /?, Pg^ ^ 7, Pg2 = 6, Pg^g^ = £, Pg^g^ = (, P g^ = ?? ' 

Noncommutative polynomials form an algebra K(zi, . . . , Zd) over K, often called 
the free associative algebra on d generators zi,. . . ,Zd- Notice that we can evaluate 
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a noncommutative polynomial p e IK(zi, . . . , Zd) on a d-tuplc Z — {Zi, . . . , Zd) of 
n X n matrices over K, for any yielding a n x n matrix p{Z). 

A non-zero polynomial can vanish on tuples of matrices of a certain size. E.g., 
p ~ Z1Z2 — Z2Z1 vanishes on pairs of 1 x 1 matrices (scalars), and 



vanishes on pairs of n x n matrices (here 5^+1 is the symmetric group on n + 1 
elements); see [JTl Proposition 1.1.37 and Exercise 1.4.11 on page 104] and [IHl 
Theorem 7]. However, if p{Z) = for all d-tuples Z of square matrices of all 
sizes, then necessarily p is the zero polynomial. More precisely, if p{Z) = for 
all d-tuples Z oi n x n matrices, then degp > 2n; this follows by applying to p 
a multilinearization process to reduce to the case of a polynomial of degree 1 in 
each indeterminate, and then evaluating the resulting multilinear polynomial on a 
staircase of matrix units, see [JTl page 6 and Lemma 1.4.3]. We notice that a much 
stronger statement appears in [1]; in particular, the intersection of the kernels of 
the matrix evaluations p{Z) for all d-tuples of n x n matrices is zero for n large 
enough compared to the degree of p. 

We next define (scalar) noncommutative rational expressions by starting with 
noncommutative polynomials and then applying successive arithmetic operations — 
addition, multiplication, and inversion. We emphasize that an expression includes 
the order in which it is composed and no two distinct expressions are identified, e.g., 
(2:1) + (— zi), (—1) + (((zi)~^)(zi)), and are different noncommutative rational 
expressions. A noncommutative rational expression r can be evaluated on a d-tuple 
Z oi n X n matrices in its domain of regularity, domr, which is defined as the set of 
all d-tuples of square matrices of all sizes such that all the inverses involved in the 
calculation of r{Z) exist. E.g., if r = {ziZ2 — Z2Zi)^^ then domr = {Z ~ {Zi, Zq) : 
det(ZiZ2 — ^2.^1) 7^ 0}. We assume that domr 7^ 0, in other words, when forming 
noncommutative rational expressions we never invert an expression that is nowhere 
invertible. 

Two noncommutative rational expressions r\ and r2 are called equivalent if 
domri ndomr2 ^ CQ and ri{Z) = r2{Z) for all d-tuples Z e domri ndomr2. E.g., 
the three different noncommuting rational expressions in the paragraph above are 
equivalent. For another example, easy matrix algebra shows that ri = z\Z2{z\Z2 — 
Z2Zi)~^ and r2 = 1 + 2:221(2:12:2 — Z2Zi)~^ are equivalent. 

We define a noncommutative rational function to be an equivalence class of non- 
commutative rational expressions. We usually denote noncommutative rational 
functions by German (Fraktur) letters. 

Notice that, unlike in the commutative case, the "minimal complexity" of a non- 
commutative rational expression defining a given noncommutative rational function 
can be arbitrarily high; there is nothing similar to a coprime fraction representation. 



This requirement is in fact superfluous. For any noncommutative rational expression r, it 
turns out 1341 Remark 2.3] tliat tlie evaluation of r on d-tuples of n X n generic matrices — see 
the discussion following Proposition 12.11 below — is defined for all sufficiently large n. In the 
case where K is an infinite field it follows that domn t is Zariski dense for all sufficiently large n; 
therefore for any two noncommutative rational expressions ri and r2, domri H domr2 7^ 0. The 
case where the field K is finite can be handled as in the proof of Proposition [2Tl] 
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It turns out that any nonzero noncommutative rational function is invertible. 
This follows from the following result which is essentially well known and is non- 
trivial already in the case where r in the statement is a noncommutative polynomial. 

Proposition 2.1. IJ r is a noncommutative rational expression and detr(Z) = 
for all Z £ domr then r is equivalent to zero. 

Proposition 12 .11 means that noncommutative rational functions form a skew field 
— a skew field of fractions of the ring of noncommutative polynomials. We remind 
the reader that a skew field, also called a division ring, is a ring with identity in 
which every nonzero element has a multiplicative inverse; it is therefore similar to 
a field, except that the multiplication is not assumed to be commutative. A skew 
field of fractions of a given ring is a skew field containing the ring and generated by 
it in the sense that no proper skew subfield contains the ring. If a noncommutative 
integral domain TZ satisfies the so called right Ore condition, 

(2.1) \/a,beTZ,b^0 3c,den,d^Q: ad = bc, 

then one can construct a skew field of fractions analogously to the commutative 
case as the ring of right quotients, i.e., of formal fractions cd~^, d ^ 0. In general, 
a skew field of fractions of a noncommutative integral domain might or might not 
exist. We refer to [H Sections 0.5-0.6 and Chapter 7], [20], [31 Chapter 4], [13 
Section 1.7 and pages 105, 107-108] for more material on skew fields of fractions 
and noncommutative localization. If a skew field of fractions exists, it might not 
be unique; we will discuss this in more detail later in this section. 

Before proceeding to the proof of Proposition 12. 1[ we notice that for a given 
matrix size n, we can view a noncommutative polynomial or a noncommutative 
rational expression in d noncommuting indeterminates as a n x n matrix-valued 
polynomial or rational function in commuting matrix entries. More formally, let 
T — (Ti, . . . , Td) be a c?-tuple of n x n matrices whose entries {Ti)jk are dn^ com- 
muting indeterminates; Ti, . . . , Td are often called generic matrices. For a noncom- 
mutative polynomial p or for a noncommutative rational expression r, we define 

Pn = p(T) e K"><" [{T,),k ■.i^l,...,d; j,k^l,...,n] 

and 

r„ = r(r) e K"^"((r,),fc : z = l, . . . , d; j, fc = l, . . . , n). 
Note that r„ is defined only for n in a subset J\fr C N such that all the inversions 
involved in the calculation of the rational matrix-valued function r{T) exist. E.g., 
if r = (ziZ2 - Z2Zi)-'^ then = N \ {1}. 

We make two useful remarks. First, if K is an infinite field then n G J\fr if and 
only if dom„ r := domr n (K"^")'' is nonempty. This may fail in the case of a 
finite field K. E.g., if K = Z/pZ and r{zi) — {z\ — zi)^^ then 1 G A/",., however 
donii r — $. The problem is, of course, in that there are nonzero polynomials over 
K which vanish identically. 

Second, we define the extended domain of regularity, edomr, of a noncom- 
mutative rational expression r as follows: for each matrix size n, edom„ r := 
edomr n (K"^")'' is the domain of regularity of the rational matrix- valued func- 
tion r„. (We set edom„ r = if n ^ Mr-) Here the domain of regularity of a 
rational matrix-valued function of commuting indeterminates is defined to be the 
intersection of the domains of regularity of its entries; the domain of regularity of a 
scalar rational function is the complement of the zero set of its denominator in the 
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coprime fraction representation. In general, domr C edomr. As a silly example, 
take r = zizj"^. Then domr = {Zi: detZi ^ 0}, however J\fr — N, r„ = /„ for 
each n, so that edomr consists of all square matrices Zi over IK of all sizes. For a 
more conceptual example, see the end of this section. 

Proof of Provosition \2.1[ We consider a subring Gn = {Pn : P G . . . , Zd)} of 

]K"^"[(Ti)jfc], which is often called the ring of generic matrices. Let £>„ be the ring 
of central quotients of Gn, i-e., the ring of formal fractions PQ~^ with P,Qg G„ 
and Q a regular central element (see [47l Section 1.7] for details). By a theorem 
of Amitsur ^71 Theorem 3.2.6], Dn is a skew field. Assuming d > 1, any central 
element Q € Gn is central in the whole ring K"^"[(rj)jfe] (see [HI Exercise 2.4.2 
on page 149]). Hence, any such nonzero Q is a scalar n x n matrix of polynomials 
in {Ti)jk; in particular, Q is invertible in K'^^^(^{Ti)jk) ■ It follows that the skew 
field Dn is a subring of K"^" ((Ti)^^) (clearly, this is also true in the case of d = 1, 
where Gn is commutative, every element is central, and every element is invertible 
in K"^"((Ti)j7£)). Therefore, rn € Dn for any noncommutative rational expression 
r with n € Nr- If det r„ = then r„ is not invertible in K"^"((Ti)jfc). On the 
other hand, since Dn is a skew field, this can happen only if r„ — 0. 

Assume now that detr(Z) = for all Z £ domr. We claim that detr„ = for 
each matrix size n S Mr', by the preceding paragraph, this will imply the conclusion 
of the proposition: r{Z) — for all Z E dom„ r, n E Mr, i.e., r is equivalent to 
zero. If the field K is infinite then the assumption det r{Z) — for all Z E dom„ r 
is simply equivalent to detr„ = (as a rational function of [Ti)jk). 

Let now K be a finite field, and assume that detr„ ^ 0. Then there exists a 

finite extension K of K, say of degree m, and Z E ^K"^"^ such that Z E domr 

(over K) and detr(Z) ^ 0. Since K is finite, it is perfect, therefore every finite 
extension is separable, and thus admits a primitive element a, i.e., K = IK(a) (see 
[391 Theorem 7.6.1 and Corollary 7.7.8]). Then I, a, a™-^ is a basis for K over 
K, and we can define a ring homomorphism (f>: K ^ j^mxm setting cj){a) — A 
where A is a m x m matrix over K whose minimal polynomial coincides with that 
of a. (j> induces a ring homomorphism = idK^x™ (8)k ((> from K"^" ^ K"^" (8)k K 
to K'^^" (g)K K"^^ = T^nmxnm^ rj,^^ ^^^^ ^^int 0„ is a homomorphism implies that 
(0„(Zi), . . . ,(/)„(Z<j)) e domr and 

ri(t)niZi), . . . ,(f>n{Zd)) = (t)n{r{Zi,...,Zd)). 

For any X — [xij] E K"^", we claim that 

det(/)(detX) = det(/)„(X), 

which boils down to 

det I ^ sign(7r)0(xi^(i)) • • • 4>ixn^(n)) 1 = det [0(%)],j=i,...,„ • 

\7rGS„ / 

Since the matrices 4>{xij) commute, the last formula follows from a well known 
identity for the determinant of a block matrix with commuting blocks; see, e.g.. 
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[37]. Finally, 

detr (^(/jniZi), . . . ,(j)n{Zd)j =det(/)„ (j{Zi, . . . , Za)^ 

= det(/)(detr(Zi,...,Zd)) . 

Since detr(Z) ^ 0, and </) is a ring honiomorphism whose domain is a field, 
(p (^detr{Z)j is invertible. Therefore, 

detr (0„(Zi),...,0„(Zrf)) =det0(detr(Z)) ^0, 

which contradicts the assumption that det7-(Z) = for each Z e donir. □ 

The proof of Proposition l2 . 1 [ implies two interesting facts. First, while the "mini- 
mal complexity" of a noncommutative rational expression defining a given noncom- 
mutative rational function can be arbitrarily high, its restriction to n x n matrices, 
for every matrix size n, is of a simple form: it is equal to the restriction of pq~^, 
where p and q are noncommutative polynomials, q being a central polynomial for 
n X n matrices. 

Second, we could have defined noncommutative rational expressions and their 
equivalence using evaluation on generic matrices or on matrices over the algebraic 
closure of K rather than evaluation on matrices over K as we did. However, we would 
obtain the same noncommutative rational expressions and the same equivalence 
relation. 

We denote the skew field of noncommutative rational functions in zi , . . . , by 
K<(;zi, . . . , Zd'^': it is often called the free skew field. 

Unlike in the commutative case, skew fields of fractions are in general not unique. 
Here is an example of infinitely many embeddings of K{zi, Z2) into skew fields (see 
[22] and [ini Exercise 7.2.10 on page 258]). Consider the polynomial ring in one 
indeterminate K[t] with the endomorphism q;„ {n — 2,3, . . .) induced by a„(i) — 
and let K[i][a;;a„] be the skew polynomial ring over 'K[t\ determined by «„. The 
elements of K[t][a;; «„] are polynomials over K in t and x with the indeterminates 
t and X satisfying the commutation relation tx = xt"'. Like any skew polynomial 
ring, K[i][x; is a right Ore ring, i.e., satisfies (|2.ip . and can be embedded into its 
skew field of right quotients, K{t){x; q;„). The elements of K{t){x; an) are rational 
functions over K in t and x with tx = xV^ . See [16l Section 0.8] for details on skew 
polynomial rings. Let y = xt. It turns out that for any noncommutative polynomial 
p in two indeterminates, p{x, y) 7^ 0. This can be verified directly by showing that 
distinct noncommutative monomials in x and y yield distinct monomials of the 
form x^t^; it can also be deduced from the result of [33], since the two left ideals in 
K[t][j:; an] generated by x and by y have a trivial intersection. We thus obtain an 
embedding p p(x, y) of K(zi,Z2) into K[i] [a;; a„] and, therefore, into K(t)(a;; a„). 
It is not hard to see that these embeddings are not isomorphic for distinct values 
of n. 

It is instructive to observe that in the embedding ]K(zi,Z2) ^ K(i)(a;;Q;„) the 
images x and y of zi and Z2 satisfy the rational identity x~^yx = y{x~^y)"'~^ (since 
tx = xt" and y — xt, hence t = x~^y). This is in contrast to the free skew field 
K<(;zi, Z2^, where zi, zi satisfy no nontrivial rational identities. 

The non-uniqueness issue leads us to the notion of the universal skew field of 
fractions. A skew field of fractions K of a ring Tl is called a universal skew field 
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of fractions if for every homomorphism : 7?. — s> L to a skew field L there exists a 
subring Kg C K containing TZ and a liomoniorphisni 9 : Kq — s> L extending such, 
that the following holds: 
(2.2) 

/or every x ^ in Kq its inverse x ^ belongs to Kq if and only if 9{x) ^ 0. 

Furthermore, the extension 9 is unique in the following local sense. Let 9' : Kq L 
be another extension satisfying (|2.2p . Then there exists a subring Kq C Kq n Kq 
containing TZ such that 9 and 9' agree on Kq, and the extension 9 : Ko L defined 
by 9 :— 9\^^_^ — 9'\-^^ satisfies (|2.2|) . It is straightforward to see that a universal 
skew field of fractions is unique (when it exists) up to a unique isomorphism. 

The following result is essentially the first fundamental theorem of Amitsur on 
rational identities. 

Proposition 2.2. ]K<(;zi, . . . , is the universal skew field of fractions of the ring 
K{zi, . . .,Zd). 

Proof. Let (/>: K{zi, . . . , Zd) L be a homomorphism to a skew field L. We notice 
that for any skew field D over K, a rational expression r can be evaluated on a 
d-tuple X = (xi, . . . , Xd) of elements of D provided all the inverses involved in the 
calculation of rix) exist, i.e., all the elements to be inverted are nonzero. Therefore, 
we can define 9{r) as the evaluation r(0(zi), . . . , (t>{zd)) whenever this is possible. 

We claim that if r is equivalent to then r(0(zi), . . . , 4){zd)) is either or unde- 
fined. Indeed, r being equivalent to means that for every matrix size n, r is either 
or undefined on the skew field of fractions Dn of the ring of generic matrices G^, 
see the proof of Proposition 12.11 In other words, r is a rational identity for Z3„, 
n = 1,2,..., hence (see [47, Theorem 8.3.3 and Corollary 8.2.16]) r is a rational 
identity for any skew field over K. In a little bit more details, the fact that r is a 
rational identity for Z?„, n — 1,2, . . . (or just for a sequence Dn , Uj — >■ oo), implies 
by a simple ultraproduct construction in the proof of Corollary 8.2.16 in j47] that r 
is a rational identity for a skew field D that is infinite dimensional over an infinite 
center; hence by the first fundamental theorem of Amitsur |4.7i Theorem 8.2.15], r 
is a rational identity for any skew field over K. At any rate, r is a rational identity 
for L, hence, r((/)(zi), . . . , 4i{zd)) is either or undefined. 

We can now define Kq to consist of all noncommutative rational functions r 
that can be represented by noncommutative rational expressions r such that 9(r) is 
defined, and we set 9{v) :— 9{r). It is clear that Kq is a subring of K<^zi, . . . , Zd'i^ 
containing ]K{zi, . . . , Zd), that 9: Kq — !■ L is a homomorphism extending cf) and 
satisfying (|2.2p . Furthermore, 9 is the only extension oicj) to Kq. Hence, if 9' : Kq — ^ 
L is another extension of cj) satisfying (|2.2p . then 9 coincides with 9' on Kq = Kq H 
Kq, and it is obvious that the extension 9: Kq — > L defined by 9 ■.=^ 9\^ ~ ^'Ik 
satisfies (g^l). ° □ 

Finally, we introduce matrix-valued noncommutative rational expressions and 
matrix-valued noncommutative rational functions. We start with matrix-valued 
noncommutative polynomials (having matrix rather than scalar coefficients) and 
use tensor substitutions for evaluations on tuples of matrices. E.g., if 

P = P0 + Pg,Z, + P,,Z2 + Pg2zf + Pg 
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is a matrix-valued noncommutativc polynomial of total degree 2 with coefficients 
Pw G KP^"*, |w| < 2, then for a d-tuple Z ~ {Zi, . . . , Zd) oi n x n matrices over K, 



P{Z) = PQ(E)In+Pg,(E)Zi+Pg,(E)Z2+Pg2(E)Zf+P, 



9291 -o-^z^l I ^ 



with ] 



pn X gn 



i.e., with 



We will often use the canonical identification of K^'^'?(X)]K 
p y. q block matrices with n x n block entries. Thus we will often view P{Z) above 
as a pn x qn matrix. We define matrix- valued noncommutativc rational expressions 
by starting with matrix- valued noncommutativc polynomials and applying succes- 
sive matrix arithmetic operations — addition, multiplication, and inversion, and 
forming block matrices: a pi x qi matrix of p2 x q2 matrix-valued noncommuta- 
tive rational expressions is a piP2 x 9i92 matrix- valued noncommutativc rational 
expression. The domain of a matrix-valued noncommutativc rational expression 
R, domi?, consists of all d-tuples Z of square matrices of all sizes such that all 
the inverses involved in the calculation of R{Z) exist. E.g., consider the 1x1 
matrix-valued rational expression 



Ri = [1 0] 



-Z2 



with 



domi?i = <^ {Zi,Z2) : det 



-Z2 


-1 


'1 









-Zi 


-Z2 


-Z2 


I - 





7^0 



Notice that a 1 x 1 matrix-valued noncommutativc rational expression is not nec- 
essarily the same as a scalar noncommutativc rational expression, since it may 
involve, as in this example, intermediate matrix operations. 

Equivalence of matrix-valued noncommutativc rational expressions, and matrix- 
valued noncommutativc rational functions as equivalence classes, are defined as 
in the scalar case. E.g., using a standard Schur complement calculation, we can 
observe that the 1x1 matrix- valued noncommutativc rational expression i?i above 
is equivalent to the following two scalar noncommutativc rational expressions. 



r2 



(1 - Zi - Z2(l - Zi) ^Z2y 



and 



-Z^\l - Z1)(Z2 - (1 - Z,)Z^\1 ~ Z,)r\ 



with 



domra - {(Zi, Z2) : det(/ - Zi) ^ 0, det(/ - Zi - ^2(1 - Zi)-^Z2) ^ 0} 



and 

domrs = {{Zi, Z2) : dct(Z2) ^ 0, det(Z2 - (/ - Zi)Z^\l - Zi)) ^ 0}. 

It is not a priori clear whether & p x q matrix-valued noncommutativc rational 
function is the same thing as a p x q matrix of (scalar) noncommutativc rational 
functions; the question is whether any pxq matrix- valued noncommutativc rational 
function can be represented hy a pxq matrix of scalar noncommutativc rational ex- 
pressions. It turns out that this is true, because noncommutativc rational functions 
form a skew field; see [Ml Remarks 2.16 and 2.11] for details. 
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We define the domain of regularity of a matrix-valued noncommutative rational 
function as the union of the domains of regularity of all matrix- valued noncom- 
mutative rational expressions representing this function, i.e., 



We emphasize that even for the case of a (scalar) noncommutative rational func- 
tion, we define its domain using all 1 x 1 matrix-valued noncommutative rational 
expressions representing the function, not just the scalar ones. E.g., in the exam- 
ples above, it is easily seen that domr2 and domrs are both properly contained 
in domi?i; so, if r is the corresponding noncommutative rational function, then 
domr D domi?!. In fact, the result on the singularities of minimal realization (to 
be discussed in Section [3]) implies that domr = domi?i. See [SH Remark 2.11] for 
additional discussion and references. 

We can also evaluate a matrix-valued noncommutative rational expression R on 
generic matrices as in the discussion preceding the proof of Proposition 12.11 and 
introduce a subset Afji C N where the evaluation is defined and the extended domain 
edom R. We then define A/th and edom *H, the extended domain of regularity of a 
matrix-valued noncommutative rational function IH, by 



We notice that while in general, for i? e 5H, domi? C dom*H, it is always the case 
that edom R — edom provided that J\fji — A/Jh ; see |34l Section 2] for additional 
discussion. 

3. Realization theory for noncommutative rational functions 

It is a bitter experience that the constellation of foundational facts underlying 
the classical Kalman realization theory for ID systems collapses for rational func- 
tions of several commuting variables and commutative multidimensional systems. 
It is all the more amazing that these facts do hold, with obvious modifications, 
in the noncommutative setting. Noncommutative systems of the form (j3.ip below 
were first studied by Ball-Vinnikov [S] in the conservative setting, in the context of 
operator model theory for row contractions (Popescu [43l EH HH US] ) and of repre- 
sentation theory of the Cuntz algebra (Bratelli-Jorgensen [13] and Davidson-Pitts 
[21]). On the other hand, noncommutative realizations very similar to p.2p were 
considered much earlier in the theory of formal languages and finite automata in the 
work of Kleene, Schiitzenberger and Fliess [3S1 UHl [21] ■ A comprehensive study of 
noncommutative realization theory appears in Ball-Groenewald-Malakorn [H [6l [5] ; 
these papers give a unified framework of structured noncommutative multidimen- 
sional linear systems for different kinds of realization formulae. We also mention 
the paper by Ball-Kaliuzhnyi-Verbovetskyi [7] where an even more general class of 
noncommutative systems (given though in a frequency domain) was described and 
the corresponding dilation theory was developed. 

A noncommutative multidimensional system is a system with evolution along the 
free semigroup J-d on d letters gi, . . . , rather than along the multidimensional 
integer lattice Z"*. An example of system equations with evolution along Td is 
given by a noncommutative Fornasini-Marchesini system (see [27] for the original 




A/'m = y Afl, edomfH= (J cdomi?. 
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commutative version): 

{x{giw) = Aix{w) + Biu{w), 

x{9dw) = Adx{w) + Bdu{w), 

y{w) = Cx{w) + Du{w), 

Applying to the system equations p.ip an appropriately defined formal non- 
commutative z-transform and under the assumption that the state of the system is 
initialized at (so that x{0) — 0), we arrive at the input-output relation 

y{z) = Tsfm(z)w(z) 

where the transfer function is given by 

(3.2) rsFM(z) + C{I^ - Aizi AdZdj-^BiZi + ■■■ + BdZd). 

Here m is the dimension of the state space K™ where vectors x{'w) live. We see 
that the transfer function is a matrix-valued noncommutative rational function 
in noncommuting indeterminates zi, . . . , which is regular at zero, i.e., zero be- 
longs to its domain of regularity (a little more precisely, the transfer function is 
the matrix- valued noncommutative rational function defined by the matrix-valued 
noncommutative rational expression p.2p '). 

The system (j3.ip is called controllable (resp., observable) if 

span ran{^"'Bj} = K", (resp., p| ker{C^"'} = {0}). 

The following facts are fundamental for the noncommutative realization theory: 

(1) Every matrix- valued noncommutative rational function which is regular at 
zero admits a state space realization (j3.2p . 

(2) An arbitrary realization (j3.2l) of a given matrix-valued noncommutative 
rational function can be reduced via an analogue of the Kalman decompo- 
sition to a controllable and observable realization. 

(3) A realization (|3.2p is controllable and observable if and only if it is mini- 
mal, i.e., it has the smallest possible state space dimension, and a minimal 
realization is unique up to a unique similarity. 

(4) A minimal realization p.2p can be constructed canonically and explicitly 
from a matrix-valued noncommutative rational function by means of the 
corresponding Hankcl operator; this ties in with the fact that the Hankel 
operator corresponding to a matrix-valued noncommutative formal power 
series has finite rank if and only if the power series represents a rational 
function (an analogue of Kronecker's Theorem). 

(5) In a minimal realization p.2p . the singularities of the transfer function 
coincide with the singularities of the resolvent; more precisely, the domain 
of regularitjQ of the transfer function (|3.2p is exactly 

{(Zi, . . . , Zrf) : det(/ - ® Zi Ad®Zd)i-^}. 

For the proofs of items (l)-(4), including missing details and exact references 
to the earlier literature, we refer to |3] where these facts are established in a more 
general setting of structured noncommutative multidimensional systems. 



^In fact, this is also the extended domain of regularity of the transfer function. 
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As for item (5), it is amazingly difficult to prove "by hands"; the usual proofs for 
d = 1 use the Hautus test for controllability / observability, but this is no longer 
available. A proof appears in [341 using noncommutative backward shifts which 
are a particular instance of the difference-differential calculus for noncommutative 
rational functions. This is a special case of the difference-differential calculus for 
general noncommutative functions, which are functions on tuples of square matrices 
of all sizes which respect direct sums and simultaneous similarities. The forthcom- 
ing basic reference is |35) . The difference-differential calculus for noncommutative 
rational functions can be developed in a more straightforward manner than in the 
general case, and we will do this later in Section |4l 

Another important example of a structured noncommutative multidimensional 
system is a noncommutative Givone-Roesser system (for the original commutative 
version of these systems, see [28]): 
(3.3) 

xiigiw) = Aiixi{w) -\ \- AidXd{w) + Biu{w), 

Xd{gdw) = Adixi{w) -\ \- AddXd + Bdu{w), 

y{w) = Cixi{w) -\ \- CdXd{w) + Du{w), 

Here xj G K™^, i.e., the state space has d components: = © ■ • ■ © K™"*. 
The transfer function of the noncommutative Givone-Roesser system is given by 

(3.4) T^Gn{z) = D + C{Im - A(z)^)-iA(z)B, 

where A is a c? x d block matrix with blocks Aij G IK™*^'™^, and A(2:) is a d x d 
block diagonal matrix, with matrix- valued noncommutative monomials Im-iZi^ . . . , 
I-md^d on the diagonal. 

The system (j3.3|) is called controllable (resp., observable) if 

span ran{P,A'"S} = K"^ , (resp., p| ker{C^'"|K-. } = {0}), j = 1, . . . , d. 
Here P, is the orthogonal projection of the state space K™ onto its j-th component 

As we already mentioned, items (l)-(4) above hold for arbitrary structured non- 
commutative multidimensional system realizations, in particular for the noncommu- 
tative Givone-Roesser realization p.3|) . The result in [34|, i.e., item (5), has been 
proved for a much more general class of realizations than p.2p : however, this class 
does not cover all structured noncommutative multidimensional system realizations, 
and we conjecture that the result might fail for noncommutative Givone-Roesser 
realizations. 

On the other hand, the symmetry appearing in Givone-Roesser system equa- 
tion!^ makes Givone-Roesser realizations more suitable for problems where this 
symmetry is essential. For general structured noncommutative multidimensional 
systems, basic arithmetic operations on transfer functions (sum, product, inver- 
sion) correspond to certain operations on systems, in the same manner as it occurs 
in the classical ID case. For noncommutative Givone-Roesser systems, we have 



■^In the case of K = C, the adjoint system S*''^ has the same form as T,^^ , but with switching 

[Aij] co\[ 



the input and output spaces and replacing the coefficient block matrix 
adjoint. 



row[Ci] D 



by its 
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that, in addition, the adjoint of the transfer function is the transfer function of the 
adjoint system. Exploiting these correspondences, one can study noncommutative 
rational functions with certain symmetries in terms of their realizations. 

In the paper by Alpay-Kaliuzhnyi-Verbovetskyi [2j , classes of matrix- valued non- 
commutative rational functions with various symmetries were studied in terms of 
their Givone-Roesser realizations. A sample result from is a version of the so- 
called lossless bounded real lemma (cf. [5 for the general bounded real lemma in 
the noncommutative setting). Let F he a q x q matrix- valued noncommutative ra- 
tional function over the field C which is regular at zero. Let J = = J* £ C'^"?. 
Then F is called matrix- J -unitary on the set J7d of d-tuples of skew-Hermitian nxn 
matrices, n — 1,2,... (which is a noncommutative analogue of the imaginary axis 
of the complex plain) if 

(3.5) F{Z)iJ ® In)FiZy = J ® In (ZeJd) 

at all points Z E J'df) dom F. Suppose that F is a q x q matrix- valued noncommu- 
tative rational function over C which is regular at zero, and let p.3|) - p.4|) be its 
minimal noncommutative Givone-Roesser system realization. Then F is matrix- J- 
unitary on J7d if and only if 

(a) D is J- unitary, i.e., DJD* — J; 

(b) there exists an invertible Hermitian solution H = diag(iJi, . . . ,Hd), with 
Hj G (£mjxmj ^ -j-j^g Lyapunov equation 

A*H + HA = -C*JC, 

and 

B = -H-^C*JD. 

This matrix H is uniquely determined by a minimal realization p.3p - p.4p . and for 
this realization it is called the associated structured Hermitian matrix. Moreover, 
F is matrix- J -inner, i.e., in addition to p.5p . F is J-contractive on the set of all 
d-tuples of n X n matrices Zj such that Zj -f Z* > 0, n = 1, 2, . . . (this set is 
a noncommutative analogue of the right half- plane), if and only if the associated 
structured Hermitian matrix is positive definite. 

4. Difference-differential calculus 

In this section we develop the difference-differential calculus for noncommutative 
rational functions and discuss various special cases and applications: directional 
derivatives, backward shifts, finite difference formulae, higher order difference- 
differential operators, and connections with formal power series. 

The difference-differential calculus for noncommutative rational functions is based 
on difference-differential operators, 

Aj : K ^ K (g) K, j = 1, . . . ,d, 

which are noncommutative counterparts of both partial finite difference and partial 
differential operators; here K = K<(:zi, . . . , Zd^- We extend A^- to matrix-valued 
noncommutative rational functions by applying these operators entrywise; we re- 
mind the reader that a matrix- valued noncommutative rational function is the same 
as a matrix of (scalar) noncommutative rational functions. We thus have 
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Our strategy will be to define Aj on matrix- valued noncommutative rational expres- 
sions recursively, starting with matrix-valued noncommutative polynomials, and 
postulating linearity and an appropriate version of the Leibniz rule. We then check 
that equivalence is preserved, thus we can define Aj on matrix-valued noncommu- 
tative rational functions. 

To define Aj on matrix- valued noncommutative rational expressions, we will 
need to introduce matrix- valued noncommutative rational expressions in two tuples 
of noncommuting indeterminates, zi, . . . , zj, and z[, . . . , z'^. They are obtained 
by applying successive matrix arithmetic operations to tensor products of matrix- 
valued noncommutative rational expressions in , . . . , zj, and in 2:^ , . . . , and 
forming block matrices. More precisely. 

Definition 4.1. (1) If i? and R' are p x q and p' x q' matrix- valued noncom- 
mutative rational expressions in zi, . . . , Zd and in z[, . . . , z'^, respectively, 
then R(S)R' is a pp' x qq' matrix- valued noncommutative rational expression 
in zi, . . . , Zd and z[, . . . , z'^, with dom{R (g) R') = dom ii x domii', and 
the evaluation is defined by 

{R ® R'){Z, Z') = R{Z) «) R'{Z'). 

Here for Z e (K"><")'^ and Z' e (^1K"'><"')'^ we have R{Z) e 1Kp><« O 1K"><", 

R'{Z') G KP'^«' (g) K"'^"' and R{Z) ® R'{Z') e K^p'^'?'?' (g) K"^" K"'^"' 
via the canonical identification of K^^? Kp'^«' with KpJ''^««'. 

(2) If i?i and i?2 are px q matrix-valued noncommutative rational expressions 
in two tuples of indeterminates, then so is Ri + R2, dom(iii -|- R2) = 
domi?i n dom_R2, and the evaluation is given by 

(i?i + R2){Z, Z') = Ri{Z, Z') + R2{Z, Z'). 

(3) If Ri and R2 are px q and qx r matrix-valued noncommutative rational 
expressions in two tuples of indeterminates, then R1R2 is p x r matrix- 
valued, dom(iiiii2) = domiii fi domi?2, and the evaluation is given by 

{RiR2){Z,Z') = Ri{Z,Z')R2{Z,Z'). 

(4) If i? is a p X p matrix-valued noncommutative rational expression in two 
tuples of indeterminates, and det R{Z, Z') does not vanish identically on 
domii, then so is R~^, 

domJ?-^ = {{Z,Z') G domR: deti?(Z,Z') ^ 0}, 

and 

R-^{Z,Z') = R{Z,Z')-^. 

(5) If -Rob, a = 1,...,P2, b = l,...,q2, arc pi x qi matrix-valued noncom- 
mutative rational expressions in two tuples of indeterminates, then R = 
[-Raf,]o^i^..._P2;6=i,....72 isPiP2 X qiq2 matrix- valucd, 

dom R= Pi dom Rat 

a=l,...,P2; fe=l,...,92 

and 

ii(Z,Z') = [ii„6(^,^')]„=l,...,p,;6=l,...,,.- 
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We notice that for a p x q matrix- valued noncommutative rational expression 

R in two tuples of indeterminates and for Z e (K"^") and Z' G ( K" ^" j with 

{Z,Z') e domi?, the evaluation e Kp^«®K"^"®K"'^"' . We will often use 

the canonical identification of K^^? K"^" ® K"'^"' = Kp^« ® (k"^" (g) K"'^"') 

with KP""'^?""' (cf. pageEl). Thus we will often view R{Z,Z') as a pnn' x qnn' 
matrix. (An alternative interpretation of the values R{Z, Z') as linear mappings 
will be considered later — see the discussion preceding Theorem l4.8n 

Two p X q matrix- valued noncommutative rational expressions, Ri and i?2, in 
two tuples of indeterminates are called equivalent if domi?i n domi?2 ^ and 
Ri{Z, Z') = R2{Z, Z') for aU pairs of d-tuples {Z, Z') in dom i?i n dom i?2 . It would 
be natural to define matrix- valued noncommutative rational functions in two tuples 
of noncommuting indeterminates as the corresponding equivalence classes. Doing 
this in a meaningful way requires analogues of Propositions 12 . II and 12.21 see [T51IT5] 
for related issues. Here we restrict ourselves to a relatively simple situation. 

Theorem 4.2. Equivalence classes ofpxq matrix-valued noncommutative rational 
expressions in two tuples of indeterminates, which are formed by using only the rules 
(1), (2), (3), and (5) in Definition \4-l\ are in a natural one-to-one correspondence 
with p X q matrices over K ® K. 

Proof. Since a matrix- valued noncommutative rational expression is equivalent to 
a matrix of scalar noncommutative rational expressions, it is clear that any matrix- 
valued noncommutative rational expression in two tuples of indeterminates which 
is formed by using only the rules (1), (2), (3), and (5), is equivalent to a matrix 
whose entries are sums of tensor products of noncommutative rational expressions. 
It only remains to show that the corresponding elements of K ® K are uniquely 
determined. Let ri, ■ ■ ■ , Vi and r'j^, . . . , be noncommutative rational functions in 
zi , . . . , Zd and in z[, . . . , z'^ represented by noncommutative rational expressions 
ri, . . . , rg and r[, . . . , , respectively. We have to show that ii ri^r[-\- ■ ■ ■ -\-rg(S)rg 
is equivalent to zero then ri + • ■ • + ® = in K K. We may assume that 
Xi, . . . , xe are linearly independent over IK, since otherwise the number of terms in 
the tensor combination can be reduced by one. 

We may assume that r[, . . . , are not all equivalent to zero, since otherwise 
there is nothing to prove. Take Z' G domr[ n ■ • • n domr^ such that r[{Z'), . . . , 
r'giZ') are not all zero. (The existence of such a Z' is established analogously to 
the reasoning in the footnote on page [31 ) This implies that the matrix elements 
(r^(Z'))y, . . . , (r^(Z'))y are not all zero for some i and j. For an arbitrary Z G 
domnn- • -ndomr^, we have r1{Z)®r\{Z')^ Vri{Z)®r'i{Z') = 0, and therefore 

(r;(Z')).,ri(Z) -!-••• + {r',{Z')\,r,{Z) - 

is a nontrivial linear dependance relation for matrices r\(Z\ . . . , ri{Z). Therefore 
ri, ■ . ■ , Xi are linearly dependent, a contradiction. □ 

We define the domain of regularity, dom CH, of a matrix over K (g) K as the 
union of the domains of regularity of all matrix-valued noncommutative rational 
expressions in two tuples of indeterminates representing ^H. 

We can also evaluate a matrix-valued noncommutative rational expression R in 
two tuples of indeterminates on generic matrices Ti, . . . , Td and T{, . . . , T^, as in 
the proof of Proposition 12. ![ and introduce a subset Mr C N where the evaluation 
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is defined and the extended domain, edomi?. We then define Af<n and edomlH, the 
extended domain of regularity of a matrix over K(8)K, by taking the union over ah 
matrix- valued noncommutative rational expressions in two tuples of indeterminates 
representing U\. 

Remark 4.3. We can also introduce matrix- valued noncommutative rational ex- 

(1) (1) 

pressions in £ tuples of noncommuting indeterminates , . . . , , j — 1, . . . , 
£, analogously to Definition 14. 1[ except that R and R' in rule (1) are now matrix- 
valued noncommutative rational expressions in t tuples and in s tuples of indeter- 
minates respectively, with t + s = £. Namely, if R and R' are p x q and p' x q' 
matrix- valued noncommutative rational expressions in , , j = 1, 

t, and in z[^\ . . . , z^"'"', j = t + 1, . . . , £, respectively, then R (g) R' is a pp' x qq' 

matrix- valued noncommutative rational expression in z^ , . . . , z)[ ^ j — 1, . . . , ^, 
with dom(i? ® R') = domi? x domi?', and the evaluation is defined by 

{R R'){Z^'\ . . . , - . . . , ZW) R'{Z^'+'\. . . , ZW). 

We notice that for a p x q matrix-valued noncommutative rational expression 
i? in ^ tuples of indeterminates and for Z'^' e (]K"j^"j)'^, j = 1, £, with 
. . . , e domi?, the evaluation 

We will often use the canonical identification of 

= KP""* (g) (K"!"""! (g> {■■■(» (K"^-iX"f-i •••)) 

with - - Thus we will often view as a pni ■ ■ ■ ni x 

qni ■ ■ ■ ne matrix. We then define the equivalence of matrix-valued noncommutative 
rational expressions in £ tuples of indeterminates and show, as in Theorem 14. 2) that 
equivalence classes olpxq matrix- valued noncommutative rational expressions in £ 
tuples of indeterminates, which are formed by using only the analogues of the rules 
(1), (2), (3), and (5) in Definition 14.11 are in a natural one-to-one correspondence 
with pxq matrices over K**^ . We can now define the domain of regularity, dom CH, of 
a matrix over K*^^ as the union of the domains of regularity of all matrix- valued 
noncommutative rational expressions in i tuples of indeterminates representing 5H. 
We can also introduce a subset Mr C N where the evaluation on generic matrices of 
a matrix- valued noncommutative rational expression R'm i tuples of indeterminates 
is defined and the extended domain, edom i?; we then define A/in and edom 91, the 
extended domain of regularity of a matrix over K**^ . 

We proceed now with the definition of difference-differential operators Aj. For 
a, p X q matrix- valued noncommutative rational expression R, Aj{R) is a p x q 
matrix- valued noncommutative rational expression in two tuples of indeterminates. 

Definition 4.4. (1) For a matrix-valued noncommutative polynomial P G 
A,{P)=J2p^ E z^^z'^ ^{K{zi,...,za)(gK{z[,...,z'a)f'"'. 

wGJ^d u,v: w—UQjV 
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(2) If Ri and R2 are px q matrix-valued noncommutative rational expressions, 
then 

A,(i?i + i?2)-A,(i?i) + A,(i?2). 

(3) If i?i is a p X g matrix- valued noncommutative rational expression and R2 
is a g X s matrix- valued noncommutative rational expression, then 

Aj{RiR2) = Aj(i?i)(l (g> R2) + (i?i (E> l)Aj(i?2). 

(4) If i? is a p X p matrix-valued noncommutative rational expression which is 
not identically singular, then 

Aj{R-^) = -{R-^ (g> l)Aj{R){l (g) R-^). 

(5) If Rab, a, — 1, . . . ,p2, b — 1, . . . , 52, are pi x qi matrix- valued noncommu- 
tative rational expressions and R = [Rab]a=i pa- 6=1 92' ^^^^^ Aj(i?) = 

[i^ji^))ab]a=l,...,p^-b=l,...,q2^ 

{A,{R)U^A,{Rab). 

It is clear that dom Aj(_R) — domi? x domi?. 

We give some examples to illustrate Definition 14.41 

Example 4.5. For a (scalar) noncommutative polynomial of total degree two in two 
indeterminates 



P 



- f3zi + JZ2 + 5zl A- tz\Z2 + C,Z2Z\ -I- r\z\, 



we have 



and 



Ai(p) = ^(1 ® 1) + (5(1 (g) zi) + (5(zi ®\) + e(l ® 4) -I- C(^2 ® 1) 



A2(p) = 7(1 1) + ^{zx ® 1) -I- C(l ® + ?7(1 » 4) + ^(2:2 ® 1). 



Example 4.6. For the three equivalent 1x1 matrix- valued rational expressions 
r2, and introduced on page [8] we have (up to trivial equivalences) 



Ai(i?i)^([l O]0 1) 



A2(i?i)'-([l 0](8)1 




Ai(r2) ~ ((1 - zi - Z2(l - z^)-'z2)-' ® 1) 

■ (1 eg) 1 + (Z2(l - zi)-^ ® 1)(1 (8) (1 - 4)"^4)) 

.(l^(l_4_4(l_^;)-i4)"i) 
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A2(r2) - ((1 ~Zi- Z2{1 - zi)-^Z2)-^ (E> 1) 

• (1 «) (1 - Z[y^z'2 + Z2{1 - Zl)-^ (g, l) 

Ai(r3) ~ (z^-^ ® 1) (l ® (z^ - (1 - zi)4"'(l - ^'i)r') 

+ {z-\l - Zl) 1) ((Z2 - (1 - 2l)22"'(l - Zl))"' » 1) 

• (l « Z^^^l - Zj) + (1 - Z1)Z2"^ «) l) 

.(i®(z^-(i-zi)zr'(i-4))-^), 

A2(r3) ~ (z-i ® 1) (l ® zr') (l (1 - z;)(z^ - (1 - zi)zr\l - z'l))-') 

+ (Z2-1(1 - Zl) 1) ((Z2 - (1 - Zi)z2-l(l - Zi))-1 1) 
• (l (g) 1 + ((1 - Z1)Z2~^ ® l) (l ® Z^"\l - Zj))) 

.(i®(z^-(i-zi)zr'(i-^i))-^). 

Example 4.7. Consider a p x q matrix-valued noncommutative rational expression 

(4.1) R = C{Irn-AiZi AdZd)-'B, 

where m G N, Aj e K"'^™ (j = l,...,d), S G K"^?, C e K^^™. (This is a 
recognizable series realization of a p x g matrix- valued noncommutative rational 
function; see [36l gSl [24].) We have 

Aj(i?) = (C® 1) ■ ((/„ - Aizi AdZd)-^ ® 1) • Ajil ® 1) 

• (1 ® (/„ - ^iz; Adz'^y^) ■ (1 B), 

or up to trivial equivalences, 

(4.2) Aj(i?) ^ (C'(/™ - A^zi AdZd)-'Aj ® l) 

• (l®(/„-Aiz; Arfz:;)-!^) 

- (C(/,„ - Aizi AdZdy^ ® 1) • (1 ® Aj(/,„ - Aiz; Adz'^r^B) . 

It is not a priori clear from Definition 14.41 that A^- preserves the equivalence 
of matrix- valued noncommutative rational expressions and can be thus defined on 
matrix-valued noncommutative rational functions. This is a consequence of the fol- 
lowing key theorem that relates the evaluation Aj{R){Z, Z') to the evaluation of R 
on d-tuples of 2 x 2 block upper triangular matrices with Zj and Z'j on block diago- 
nals. We win identify K"^"«)K"'^"' with hom(K"^"', K"^"'), so that J2iAi®B, 
corresponds to the linear mapping H n> J^iAiHBi. This correspondence extends 

naturally to pxg matrices: we identify ]KP^«®IK"^"(g)K"'^"' ^ (K"^" ® IK"'^"^ 
with homflK"^"',(]K"^"y''), so that A^^'"'' ® ^P''''^ 



sponds to the linear mapping H i~> 



corre- 

j = l,...,p; k=l,...,q 

. We will 

j=l,...,p; k=l,...,q 



use the pn x pn permutation matrices P{p,n) — [i?j^]i=i,...,p; j=i,...,n where each 
Eij £ W^" has entry 1 in position and all other entries are zero. These 
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matrices allow us to change the order of factors m tensor products: A (g) B — 
P{n,p){B ® A)P(m,(j)"r for any A e IK"^" and B e K^^?. See [32l pages 259- 
261]. We also use the notation dom„ R domi? n (K"^")''. 



Theorem 4.8. Let R be a px q matrix-valued noncommutative rational expression. 

d 

Then 



Let Z e dom„ R and Z' G dom„/ R, and let W = (VFi, . . . , Wd) G (k""""') 

e dom„+„/ R 



'z w 




'Zi wi 






Wd 




Z' 




Z[ 







ZA 


) 



and 

(4.3) P{n + n',p)R 



Z W 
Z' 



P(n + ri . q) 



P{n,p)R{Z)P{n,q)^ P{n,p) ( ^ A,(i?)(Z, Z')(Vl^,) 1 P{n\q)^ 



P{n',p)R{Z')P{n',q)'^ 
Proof. We establish (|4.3I) recursively; the reasoning will also imply that 



Z W 
Z' 



dom„+„/ R. 

We first verify that (j4.3p holds when i? is a matrix- valued noncommutative poly- 
nomial P = Y^wdTa PwZ"". 



P{n + n',p)P 



Z W' 
Z' 



P{n + n',qV = ^ 



Z w 

Z' 



E 

te'i+,ii....,iie{i,....d} 












Z'n. 





z^, w,, 
K. 



■Pn. 



E 



k=l 







z'h ■ • • z^ 



p, 



£, zi k—l 

E 



, 11 , . . . ,1 



E^'^csp^, EE E z"w,z'-®p^ 

w w j—1 u,v: w—UQjV 

E ^'^ » 



P{n,p)P{Z)P{n,qy P{n,p) ^^^,{P){Z, Z'){Wj)j P{n',q)^ 

P{n',p)P{Z')P{n',qy 

If (|4.3p is true for Ri and for i?2 then it is clearly true for R = Ri + R2. 
Now, assume that (|4.3p is true for a pxq matrix- valued noncommutative rational 
expression i?i and for a q x r matrix-valued noncommutative rational expression 
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R2. Then 

P{n + n',p){RiR2 
= P{nW,p)Ri 



Z W 
Z' 



Z W 
Z' 



P{n + n'.ry 
P{n+n' ,q)^ P{n+n' ,q)R2 



Z W 
Z' 



P{n+n\ r) 



Pin,p)R,{Z)P{n,q)^ P{n,p)\^^A,iR,){Z,Z'){W,)\ P{n',q)^ 







P{n',p)Ri{Z')P{n',qy 



P{n,q)R2{Z)P{n,ry P{n,q) \ Y.^/^,{R2){Z, Z'){W,)\ P{n' ^r^ 





P{7i,p)Ri{Z)R2{Z)P{n,r)^ 




Pin',q)R2{Z')P{n',r 



Pin,p) iRiiZ)A,iR2)iZ,Z')iW,) + A,iR,)iZ,Z')iW,)R2iZ')) Pin' , r) ^ 
P{n',p)Ri{Z')R2{Z')P{n\r)^ 



P{n,p){R^R2){Z)P{n,ry P{n,p) \ Y.^A,{R^R2){Z, Z'){W,) \ P{n',r)^ 







Pin',p){R,R2){Z')P{n',r) 



Thus, (|4.3p is true for the p x r matrix- valued noncommutative rational expression 

Next, assume that (|4.3p is true for apxp matrix- valued noncommutative rational 
expression R which is not identically singular. Then 



P{n + n',p)R-^ 



Z W 
Z' 



^ P{n + n\p)'^ 



P{n + n',p)R 



Z W 
Z' 



P{n + Ji ,p)'^ 



P{n,p)RiZ)Pin,p)^ P{n,p) I X; A, (i?)(Z, I P{n',p)^ 





P{n,p)R{Z)-^P{n,p)'^ 




P{n\p)R{Z')P{n' ,p)-^ 



-P{n,p)R{Z)-^ ^^A,{R){Z,Z'){W,)j R{Z')-' P{n' ,p) 
P{n',p)R{Zr'P{n',pV 



P{n,p)R-\Z)P{n,py Pin,p) \^^^A,{R-^)iZ, Z')iW,) j P{n' .p^ 
P{n',p)R-^{Z')P{n',p)^ 
Thus, ()4.3p is true for the matrix- valued noncommutative rational expression R^^. 
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Finally, we will show that if ()4.3p is true for Rab , a = 1 , . . . , p2 , b = 1, . . . , q2, 
then it is also true for R = [Rab]- Clearly, it suffices to prove ()4.3|) for the case 
where only one block Rab is nonzero, i.e., R = Eab ® Rab, with some a, b and Eab a 
P2 X (72 matrix with 1 at the (a, b) position and elsewhere. Since {Eab(^ Rab){X) = 
{Eab<^Ipi)-{Iq2'®Rab)iX), this boils down to proving (|4.3I) for I q^® Rab- Simplifying 
the notation, we can state the problem as follows: show that if a p x g matrix- valued 
noncommutative rational expression R satisfies (|4.3p . then so does the sp x sq 
matrix- valued noncommutative rational expression Is® R for every s G N. We first 
observe that for any m, s S N and any X S dom„i R one has X ® 1^ & dom„js R 
and R{X) ® Is = R{X ® Is) — this follows directly from the definition of matrix- 
valued noncommutative rational expressions. Second, using this observation and 
the identities of the form 

P(m, sp) = P(ms,p)P{'mp, s) 
(see [3H Problem 20, page 266]), we obtain that 

P(m, sp){Is ® R){X)P{m, sq)'^ 

= P(nis,p)P{mp,s){Is ® R(X))P{Tnq,s)~'' P{ms,q)'^ 

= P{ms,p)R{X ® Is)P{ms,q)'^ . 



Then we have 

P{n + n,sp){Is®R) 



W 
Z' 



P{n + n , sq)'^ 



= P{{n + n')s,p)R 

'P{ns,p)R{Z ® Is)P{ns, q)~^ 




Z®Is 





Z'i 



P((7i + n')s,g)^ 



Pins,p) (Ejti Aj(i?)(^ ® Is, Z' ® Is){W, ® Is)) Pin's, g)^ 



P{n's,p)R{Z' i 
'P{n,sp){Is ® R)(Z)P(n,sq)^ 



'Is)P{n's,q) 







P{n, sp) 



P{n',sp){Is 



)R){Z, Z'){Wj)) P{n',sq)-^ 
'R){Z')P{n',sq)'^ 



i.e 



., Is® R satisfies (|43)) . 
The proof is complete. 



□ 



Remark 4.9. The equality (|4.3|) is the accurate statement of the fact that matrix- 
valued noncommutative rational expressions respect a block triangular matrix struc- 
ture, in particular direct sums — compare [331 formulae (2.3), (2.5), (2.12), (2.13), 
(2.19)], where the permutation matrices are missing. (This does not affect any 
subsequent arguments there.) 

Corollary 4.10. For any matrix-valued noncommutative rational expression R, 
edomAj(i?) 3 edomi? x edomi?. 

The proof is obtained by substituting generic matrices Ti , . . . , Td and T[, . . . , 
into - the details are similar to [34l Corollary 2.12]. 
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Corollary 4.11. If Ri and R2 are two equivalent matrix-valued rational expres- 
sions then Aj(Ri) and Aj(i?2) are also equivalent. 

The proof is immediate from (|4.3|) . 

Corollary 14.111 and Theorem 14.21 allow us to define AjDl for a matrix- valued 
noncommutative rational function U\: AjfH is the matrix over Kcg)K corresponding 
to the equivalence class of Aj (R) for any i? G $H. We have 

dom Aj5H 3 dom*K x dom9^, edomA^^H 3 edom?l x edomfH, 

where the second inclusion follows from Corollarv l4.10l 

Remark 4.12. We emphasize that while we always have equality domAj(i?) = 
dom(i?) X dom(i?) for a matrix-valued noncommutative rational expression R, we 
may have strict inclusion dom Aj9i 3 dom?l x dom5K. As an example, let r = z^^ . 
Then we have 

Aar = -(zfi (g) 1) • 0(1 (g) 1) • (1 Cg) zf^) - 0, 
and for the corresponding noncommutative rational function r we have 

oc 

domA2r= Y[ 

n=l 

while 

domr = {(Zi,Z2): dot Zi 7^0}. 
This example also shows that we have a strict inclusion 

edom(A2r) 3 edom(r) x edom(r). 

We proceed to describe some of the many facets of the noncommutative difference- 
differential operators. 

4.1. Directional derivatives. Evaluating Aj{R), j — 1, . . . ,d, at Z' = Z yields 
the differential or the directional derivatives of a matrix-valued noncommutative 
rational expression R at Z: 

d , 

(4.4) ^ A,(i?)(Z,Z)(W^,) = -R{Z + tW)l^^. 

Here Z G domi? and we view R{Z + tW) as a matrix- valued rational function in 
the indeterminate t. The formula (|4.4|) follows easily from the recursive definition 
of Aj and the corresponding properties of the derivative of a matrix- valued rational 
function in one indeterminate. Of course, one can also write the analogue of (14. 4p 
for matrix-valued noncommutative rational functions. 

4.2. Backward shifts. Let i? be a matrix-valued noncommutative rational ex- 
pression which is regular at zero, i.e., G domi R. Then 

(4.5) Aj(i?)(Z,0) = 7^Ji?(Z), Aj{R){0,Z) ^ CjR{Z), 

where TZj and Cj are the right and left backward shift operators introduced in 
[34] . The formula (|4.5|) follows easily from the recursive definitions of Aj and 
of the backward shifts. Of course, one can also write the analogue of (|4.5p for 
matrix-valued noncommutative rational functions. As an illustration of the action 
of backward shifts, for the polynomial p of Example 14.51 we have 

TZip — /3 + Szi + (z2, C-ip — (i + 5zi + ez2, 
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'R2P = ^ + ezi + riZ2, C2P = 1 + C,zi + riZ2, 
and for the recognizable series realization (|4.ip of Example 14.71 '^^ have 
TZjR^Cil^^Aizi AdZdY^AjB, C,R ^ CAj{I,n-Aizi AdZd)-^B. 

4.3. Finite difFerence formulae. For a. p x q matrix-valued noncommutative ra- 
tional expression R, and for , Z G dom„ R, we have the noncommutative finite 
difference formula 

RiZ) - R = J2 ^jiR) {Zj - ) , 

which can also be proved recursively, and extends naturally to matrix- valued non- 
commutative rational functions. 

4.4. Higher order difference- differential operators. We can iterate the differ- 
ence-differential operators Aj: we define a linear mapping 

Aj : K®^ ^ K^i^+i) 
by its action on pure tensors as 

Aj (ri (g) ■ • • r^) = ri (g) • • • ig) tg-i ® Ajxe. 
It is easy to check that Aj satisfies the following version of the Leibniz rule: 

(4.6) Aj{xx') = Aj-r • av' -f (r ® f) • A^-r', 

for all r, r' G K*^^, where the linear mapping te : K®^ K®^+^ (actually, a homo- 
morphism of K-algebras) is defined on pure tensors by 

i.i{xi (g) • • • g) r^-i (g) r^) = ri g) • • • g) r^-i (g 1 g) r^. 

Applying (|4.6p to both sides of the identity rr^^ = f , wc obtain that 

(4.7) A,(r-i) = -(r-i®f).A,r..,(r-i). 

We extend Aj and li entry wise to matrices. 

We can now define, for a word w = gig ■ ■ ■ gt^ of length £, the corresponding 
higher-order difference-differential operators 

A"' A,, ■ • • Ai, : K^^' ^ (K®(^+l))^''^ 

Example 4.13. For a matrix-valued noncommutative rational function *H defined 
by a matrix-valued noncommutative polynomial P = X^ijGJ^d ^f-^" ^'^'^ '^'-"^ ^ word 
w = gig ■ ■ ■ gi^ of length £ at most the total degree of P, we have 

A^'^n = XI (2^^^)"' g) (z(2))«= (g • • • g) 

More precisely, every term in the second sum on the right-hand side is a tensor 
product of £ -f 1 noncommutative rational functions defined by the corresponding 
noncommutative monomials; alternatively, the right-hand side (with some nesting 
of parentheses) is a matrix- valued noncommutative rational expression in ^ tuples 
of indeterminates defining A'^'lH. 
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Example 4.14. For the matrix- valued noncommutative rational function defined 
by a recognizable series realization of the form (|4.1|) . we have by iterating the first 
equality in (j4.2p that for every w = gi^ ■ ■ ■ gi-^, 

A^9\ = I C{I„r - AizJ^'> AdZ^'VM,, ® 1 • • • ® 1 



• Yl 1 • • ■ 18) 1 ^ (/m - Aiz[^^ Adz'-J^y^A,^ (8) 1 (8) ■ ■ • (g 1 

\ J — 1 times + l timcsV 

• f (Im - Al^f"-'^ Arfz('+^')-li? 

\ ^ times 

We identify K"!"""! ■ • • 8) IK"«+i''"'!+i with ^-hnear mappings 



iW^...^ .(^+1)^ ..... 

/l(2)_._ jWrr. 4(^+1) 



so that ( ® ' ■ ' ) corresponds to the ^-Unear mapping 



iH,,...,H^^j:^^H,Ar^...Af^H,Ar 

i 

This correspondence extends naturally to matrices. It follows that for a p x q 
matrix- valued noncommutative rational function *H and for Z^^^^ G K^^^"!, 
^{i+i) g j^n^+ixn^+i appropriate domains, we have that 

A^«K . . . , e KP^9 C8 K"i^"i (8 • • • 8) IK"'+iX"f+i 

^ (K"i (8 • • • K"*+i x"*+i)^''', 
and for {Hi, . . . , He) (E K"ixn2 x • • • x K'^^x'^'i+i, we have that 
A'"5h(z«,...,z(^+i)) {Hi,..., He) E 



In particular, we have 

d d ,2 

5] ^ A,A,1H(Z, Z, Z){W,, W,) = ^5n(Z -I- tW) 



t=0 



i=l ]=1 

This is exactly the Hessian of d\ which plays a central role in the study of noncom- 
mutative convexity; see, e.g., [29l |3T1 [30] . Hessians and directional derivatives of 
matrix-valued noncommutative rational expressions were implemented in NCAlge- 
bra to produce a convexity checking algorithm [141 . 

Remark 4.15. We can also define the difference-differential operators on matrices 
over tensor powers of K at the level of matrix-valued noncommutative rational 
expressions in several tuples of indeterminates. For a p x q matrix- valued noncom- 
mutative rational expression i? in £ tuples of indeterminates (see Remark 14.31) , we 
define Aj{R), a p x q matrix- valued noncommutative rational expression in £ -I- 1 
tuples of indeterminates, analogously to Definition 14.41 except that in rule (1) we 
consider "pure tensors" instead of polynomials, and we modify rules (3) and (4), 
cf. ([Ml and ([HZ]). Namely, 
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(1) If i? and R' arc matrix- valued noncommutative rational expressions in t 
tuples and in s tuples of indeterminates respectively, with t + s = we set 

Aj{R(g)R') = Aj(i?'). 

(3) If i?i and i?2 ai'e matrix-valued noncommutative rational expressions of 
compatible sizes in £ tuples of intederminates, then 

A,(i?ii?2) = A,(i?i)i,(i?2) + ® l)Aj(i?2). 

(4) If i? is a square matrix-valued noncommutative rational expression in £ 
tuples of indeterminates, which is not identically singular, then 

Aj{R-^) = -{R-^ ® l)Aj{R)Li{R-^). 

Here, for a matrix-valued noncommutative rational expression i? in ^ tuples of 
intederminates, ii{R) is a matrix- valued noncommutative rational expression of 
the same size m £+1 tuples of intederminates defined by ti (R) = 1(E) R and by the 
recursive relations 

(1) H+s{R^R') = R®is{R'); 

(2) ie{Ri+R2) = Li{Ri) + Li{R2); 

(3) /,,(i?ii?2) - t,(i?i)tf(i?2); 

(4) = 

(5) Le{[Rab]) = MRab)] 

for ^ > 1. Notice that for d-tuples of matrices e (K"^ ^"^ j = 1, . . . , ^-|- 1, the 
value is the image of the value R{Z^^\ . . . , Z^'^-'^\ Z^'^+'^^) 

under the linear mapping 

K"i^"i ® ...(g)K"'-i><"*-i (g)K"*+i^"'^+i 

— ^K"!"""! ® •••(8)IK"«-i''"'-i (K)K"*''"* (g)K"*+i''"*+i, 
A^^^ ® • • • (g) yl(^) ^ «) A(^) 

extended naturally to the mapping of matrices 



We then have an analogue of Theorem 14.81 as follows. Let R he a. p x q matrix- 
valued noncommutative rational expression in £ tuples of indeterminates, let Z^^-* S 

(W^x^j f^ j =.1, £-1, Z e (K"^")'', Z' e (^K"'^"') , so that 

edomi?, 

and let We (k"^"')"*. Then 



G domi? 
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and 
(4.8) 

P{n + n',pni---ni^i)R ( Z'-^\ . 



z w 

Z' 



^ P{n + n', qni ■ ■ ■ ne^i)~^ 



P{n,pni ■ ■ ■ ne^i)R{Z(^\ Z<^'-^\ Z)P{n, qm--- 





P{n,pni ■ ■ ■ n,„i) ^ A, , . . . , Z('-^\ Z, Z'){W,)P{n', qm ■ ■ ■ n,_i)^ 

i=i 

P{n',pni ■ ■ ■ . . . , Z('-^\Z')P{n', qm ■ ■ ■ 

Here 

A,(i?)(zW, 



pxq 



and 

A,(i?)(Z(i), 



pxg 



cf. the discussions preceding Theorem 14.81 and following Example 14.141 The proof 
is analogous to the proof of Theorem 14.81 except that instead of establishing ()4.3p 
for polynomials we have to establish (j4.8l) for pure tensors. Namely, we have to 
show that if R and R' satisfy (|4.8p then so does RiS$ R'. This can be achieved using 
the identity 



= fl^R' 



z w 

Z' 



It follows from (j4.8p that Aj preserves the equivalence of matrix- valued noncom- 
mutative rational expressions in £ tuples of indeterminates and can be thus defined 
on matrices over K*^^. It also follows that if 

eedomi?, 



e edomAj(i?) 



then 



Remark 4.16. As a special case of the previous remark, we see that for a matrix- 
valued noncommutative rational function and for a word w of length A^9\ 
is the matrix over K®^ corresponding to the equivalence class of A'"(i?) for any 
R G We further conclude that 

domA^-^H D (dom5n)^+\ edomA^'D^ 2 (edom$H)^+^ 

The second inclusion follows from the last statement of Remark 14.151 The first 
inclusion follows from the equality domA^i? = (domi?)^+^ for a matrix- valued 
rational expression R, which can be proved recursively. We leave the details to 
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the reader, noticing only that for the product we have the fohowing higher order 
Leibniz rule: 

(4.9) A"'(i?ii?2)= J2 (A^li^O^l^-'-^lj-Cl®---®!® A"(i?2)). 

4.5. Formal power series. A matrix-valued noncommutative rational expression 
which is regular at zero determines a noncommutative formal power series with 
matrix coefficients. This correspondence is defined recursively by inverting formal 
power series with invertible constant term (the coefficient for z*^); see, e.g., |12| . 
Furthermore, i?i and R2 are equivalent if and only if the corresponding formal 
power series coincide, so that the noncommutative formal power series expansion 
of a matrix-valued noncommutative rational function which is regular at zero is 
well defined; see [34l Remark 2.14]. If J^wGJ^a^^^^ ^^'^ formal power series 
expansion of D\ then the formal power series expansion of Aj9l is given by 

For a proof, we use the recognizable series realization and represent ?l by a matrix- 
valued noncommutative rational expression of the form (j4.ip . Therefore, IH^ = 
CA^B. On the other hand, it follows from Example 14.71 that the formal power 
series expansion of Aj*H is given by 

So, we see that the coefficient for z" (g) z'" is exactly ^ugjv 

A similar argument using Example 14.141 shows that the formal power series ex- 
pansion of A'^IH for w = gi^ ■ ■ ■ gi-^ is given by 

E «n,,.,„.,,,...,,,„.,,(z(^))"^ ® (z^'^r ® • • • ® (z(^+i))«^+^ 

In particular, looking at the constant term of this expansion (i.e., for ui 

Ui^i = 0), we see that 

A'"$R(0,...,0)=^H„x. 



5. Conclusions 

Rational functions in noncommuting indeterminates occur in many areas of sys- 
tem theory: most control problems involve rational expressions in matrix param- 
eters. In this paper we surveyed some aspects of the theory of noncommutative 
rational functions, and provided some pointers to a rapidly growing literature. We 
discussed in some details a construction of the skew field of noncommutative ra- 
tional functions based on noncommutative rational expressions and their matrix 
evaluations. We explained its role as the universal field of fractions of the ring of 
noncommutative polynomials. We gave an outline of a noncommutative realization 
theory. Finally, we developed a difference-differential calculus for noncommutative 
rational functions. 



NONCOMMUTATIVE RATIONAL FUNCTIONS 



27 



Aknowledgements 

It is our pleasure to thank two anonymous referees and Amnon Yekutieli for use- 
ful comments and suggestions. We are also grateful to Igor Klep for the l^jKmacros 
used in the notation for the free skew field. 

References 

[1] D. Alpay and D. S. Kalyuzhnyi-Verbovetzkii, On the intersection of null spaces for matrix 

substitutions in a non-commutative rational formal power series. C. R. Math. Acad. Sci. 

Paris 339(8):533-538, 2004. 
[2] D. Alpay and D. S. Kalyuzhnyi-Verbovetzkii. Matrix- J- unitary non-commutative rational 

formal power series. In The state space method generalizations and applications, volume 161 

of Oper. Theory Adv. Appl., pages 49-113. Birkhauser, Basel, 2006. 
[3] S. A. Amitsur. Rational identities and applications to algebra and geometry. J. Algebra, 

3:304-359, 1966. 

[4] J. A. Ball, G. Grocncwald, and T. Malakorn. Structured noncommutativc multidimensional 
linear systems. SIAM J. Control Optim., 44(4): 1474-1528 (electronic), 2005. 

[5] J. A. Ball, G. Groenewald, and T. Malakorn. Bounded real lemma for structured noncommu- 
tative multidimensional linear systems and robust control. Multidimens. Syst. Signal Process., 
17(2-3): 119-150, 2006. 

[6] J. A. Ball, G. Groenewald, and T. Malakorn. Conservative structured noncommutative mul- 
tidimensional linear systems. In The state space method generalizations and applications, 
volume 161 of Oper. Theory Adv. Appl, pages 179—223. Birkhauser, Basel, 2006. 
[7] J. A. Ball and D. S. Kaliuzhnyi-Vcrbovctskyi. Conservative dilations of dissipativc multidi- 
mensional systems: The commutative and non-commutative settings. Multidim. Syst. Signal 
Processing, 19:79-122, 2008. 
[8] J. A. Ball and V. Vinnikov, Lax-Phillips scattering and conservative linear systems: a Cuntz- 

algebra multidimensional setting. Memoirs Amer. Math. Soc. 178, no. 837, 2005. 
[9] C. Beck. On formal power series representations for uncertain systems. IEEE Trans. Automat. 
Control, 46(2):314-319, 2001. 

[10] C. L. Beck, J. Doyle, and K. Glover. Model reduction of multidimensional and uncertain 
systems. IEEE Trans. Automat. Control, 41(10):1466-1477, 1996. 

[11] G. M. Bergman. Skew fields of noncommutative rational functions, after Amitsur. In 
Seminaire Schiitzenberger -Lentin-Nivat, Annee 1969/70, No. 16. Paris, 1970. 

[12] J. Berstel and C. Reutenauer. Rational series and their languages, volume 12 of EATCS 
Monographs on Theoretical Computer Science. Springer- Verlag, Berlin, 1988. 

[13] O. Bratelli and P. E. T. Jorgensen, Iterated function systems and permutation representations 
of the Cuntz algebra, Memoirs Amer. Math. Soc. no. 139, 1999. 

[14] J. F. Camino, J. W. Helton, R. E. Skelton, and J. Ye, Matrix inequalities: A Symbolic 
Procedure to Determine Convexity Automatically, Integral Equations and Operator Theory, 
46(4):399-454, 2003. 

[15] P. M. Cohn. The embedding of firs in skew fields. Proc. London Math. Soc, 23:19.3-213, 1971. 

[16] P. M. Cohn. Free rings and their relations. Academic Press, London, 1971. London Mathe- 
matical Society Monographs, No. 2. 

[17] P. M. Cohn. Universal skew fields of fractions. Symposia Math., 8:135-148, 1972. 

[18] P. M. Cohn. The universal skew field of fractions of a tensor product of free rings. Colloq. 
Math. 72(l):l-8, 1997. 

[19] P. M. Cohn. Correction to: "The universal skew field of fractions of a tensor product of free 
rings" [Colloq. Math. 72 (1997), no. 1, 1-8]. Colloq. Math. 76(2):319, 1998. 

[20] P. M. Cohn. Free ideal rings and localization in general rings. Cambridge University Press, 
Cambridge, 2006. New Mathematical Monographs 3. 

[21] K. R. Davidson and D. R. Pitts, Invariant subspaces and hyper-refiexivity for the free semi- 
group algebras, Proc. London Math. Soc. 78:401-430, 1999. 

[22] J.L. Fisher, Embedding free algebras in skew fields. Proc. Amer. Math. Soc. 30:453—458, 
1971. 

[23] M. Fliess. Sur le plongement de I'algebre des series rationnelles non commutatives dans un 
corps gauche. C. R. Acad. Sci. Paris, Ser. A, 271:926-927, 1970. 



28 D. S. KALIUZHNYI-VERBOVETSKYI AND V. VINNIKOV 

[24] M. Flicss. Matrices de Hankcl. J. Math. Pures AppL, 53(9);197-222, 1974. 

[25] M. Flicss. Sur divers produits dc scries formcUes. Bull. Soc. Math. France, 102:181-191, 1974. 

[26] E. Formanek. The polynomial identities and invariants of n x n matrices. CBMS Regional 

Conference Series in Mathematics, 78. The American Mathematical Society, Providence, RI, 

1991. 

[27] E. Fornasini and G. Marchesini. Doubly-indexed dynamical systems: state-spax;e models and 
structural properties. Math. Systems Theory, 12(l):59-72, 1978/79. 

[28] D. D. Givone and R. P. Roesser, Multidimensional linear iterative circuits-general properties, 
IEEE Trans. Computers 21: 1067-1073, 1972. 

[29] J. W. Helton. Manipulating matrix inequalities automatically. In Mathematical systems the- 
ory in biology, communications, computation, and finance (Notre Dame, IN, 2002), volume 
134 of IMA Vol. Math. AppL, pages 237-256. Springer, New York, 2003. 

[30] J. W. Helton, S. A. McCuUough, M. Putinar, and V. Vinnikov. Convex Matrix Inequalities 
versus Linear Matrix Inequalities, IEEE Trans. Ant. Control, 54(5):952-964, 2009. 

[31] J. W. Helton, S. A. McCuUough, and V. Vinnikov. Noncommutative convexity arises from 
linear matrix inequalities. J. Funct. Anal., 240(1):105-191, 2006. 

[32] R. A. Horn and C. R. Johnson. Topics in matrix analysis. Corrected reprint of the 1991 
original. Cambridge University Press, Cambridge, 1994. 

[33] A. V. Jategaonkar. Ore domains and free algebras. Bull. London Math. Soc, 1:45—46, 1969. 

[34] D. S. Kaliuzhnyi-Verbovetskyi and V. Vinnikov. Singularities of Noncommutative Rational 
Functions and Minimal Factorizations. Lin. Alg. AppL 430:869-889, 2009. 

[35] D. S. Kaliuzhnyi-Verbovetskyi and V. Vinnikov. Foundations of noncommutative function 
theory. In preparation. 

[36] S. C. Kleene. Representation of events in nerve nets and finite automata. In Automata studies. 
Annals of mathematics studies, no. 34, pages 3-41. Princeton University Press, Princeton, N. 
J., 1956. 

[37] I. Kovacs, D. S. Silver, and S. G. Williams. Determinants of commuting-block matrices. Amer. 

Math. Monthly 106(10):950-952, 1999. 
[38] J. Lambek. Lectures on rings and modules. With an appendix by Ian G. Connell. Blaisdell 

Publishing Co. Ginn and Co., Waltliam, Mass. -Toronto, Ont. -London, 1966. 
[39] S. Lang. Algebra. 2nd Edition. Addison-Wesley, 1984. 

[40] W.-M. Lu, K. Zhou, and J. C. Doyle. Stabilization of uncertain linear systems: an LET 
approach. IEEE Trans. Automat. Control, 41(l):50-65, 1996. 

[41] A. Nemirovskii. Advances in convex optimization: conic programming. Plenary Lecture, In- 
ternational Congress of Mathematicians (ICM), Madrid, Spain, 2006. 

[42] Yu. Nesterov and A. Nemirovskii. Interior-point polynomial algorithms in convex program- 
ming, volume 13 of SIAM Studies in Applied Mathematics. Society for Industrial and Applied 
Mathematics (SIAM), Philadelphia, PA, 1994. 

[43] G. Popescu, Models for infinite sequences of noncommuting operators, Acta Sci. Math. 53:355- 
368, 1989. 

[44] G. Popescu, Characteristic functions for infinite sequences of noncommuting operators, J. 
Operator Theory 22:51-71, 1989. 

[45] G. Popescu, Isometric dilations for infinite sequences of noncommuting operators. Trans. 
Amer. Math. Soc. 316:523-536, 1989. 

[46] G. Popescu, Multi-analytic operators on Fock spa<;es. Math. Ann. 303:31-46, 1995. 

[47] L. H. Rowen. Polynomial identities in ring theory, volume 84 of Pure and Applied Mathe- 
matics. Academic Press Inc. [Harcourt Brace Jovanovicli PuUishers] , New York, 1980. 

[48] M. P. Schiitzenberger. On the definition of a family of automata. Information and Control, 
4:245-270, 1961. 

[49] M. P. Schiitzenberger. Certain elementary families of automata. In Proc. Sympos. Math. 

Theory of Automata (New York, 1962), pages 139—153. Polytechnic Press of Polytechnic 

Inst, of Brooklyn, Brooklyn, New York, 1963. 
[50] R. E. Skelton, T. IwasaJci, and K. M. Grigoriadis. A Unified Algebraic Approach to Linear 

Control Design. Taylor & Francis, 1997. 



NONCOMMUTATIVE RATIONAL FUNCTIONS 



29 



Department of Mathematics, Drexel University, 3141 Chestnut Str., Philadelphia, 
PA, 19101 

E-mail address: dmitryk@inath.drexel.edu 

Department of Mathematics, Ben-Gurion University of the Negev, Beer-Sheva, Is- 
rael, 84105 

E-mail address: viniiikovOmath.bgu.ac.il 



